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TR 要 :全 面 回顾 了 非 连续 变形 分 析 (DDA ) 方法 与 数值 流 形 方法 (NMM) 在 岩石 力学 与 工程 问题 模 
〇 拟 领 域 的 研究 与 应 用 进展 。 首 先 论述 了 相关 理论 与 算法 的 研究 现状 ,包括 DDA 中 旋转 虚假 体积 脱 
胀 及 能 量 损耗 问题 的 解决 ,开裂 破坏 模拟 的 子 块 体 单 元 DDA 方法 和 节点 DDA 方法 ,以 及 颗粒 DDA 
方法 的 发 展 等 ;也 包括 对 NMM 开裂 破坏 模拟 能 力 及 双重 窗 盖 系统 的 研究 ,无 网 格 NMM 及 粒子 
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还 对 接触 相关 问题 的 研究 进行 了 专门 论述 。 随 后 ,对 DDA 与 NMM 在 岩 体 边 坡 失 稳 与 落石 地 下 将 
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表 性 研究 与 应 用 工作 进行 了 详细 论述 。 在 此 基础 上 ,凝练 和 论述 了 DDA 和 NMM 研究 与 应 用 面临 
的 主要 问题 与 挑战 ,包括 计算 控制 参数 .计算 效率 、 多 物理 场 、 三 维 方法 程序 的 商业 化 开发 等 ,并 给 
出 了 相应 的 对 策 建 议 。 
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Research progress in DDA and NMM simulations of rock 
mechanics and engineering problems 
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(1. School of Mechatronic Engineering, Southwest Petroleum University ,610500 Chengdu, China;2. Research 


Institute of Engineering Safety Assessment and Protection, Southwest Petroleum University , 610500 Chengdu, China) 


Abstract;This paper fully reviews the progress in the research and application of discontinuous deforma- 
tion analysis (DDA) and numerical manifold method ( NMM) in the field of rock mechanics and engi- 
neering problem simulations. First, the research status of related theories and algorithms is expounded, in- 
cluding the solution to the problem of spurious rotational volume expansion and energy loss in DDA, the 
sub-block element DDA method and the nodal-based DDA method for fracturing failure simulation, and 


the development of particle DDA method. Research on fracturing simulation capability and dual cover sys- 


"tems of NMM, as well as development of the meshless NMM and particle NMM are also included. Re- 


“search on simplex integration, high-order method, material constitutive model, and explicit and parallel 


"computing are also involved. Moreover, research on contact-related problems is specially discussed. Subse- 


quently, this paper expounds on the representative research and application works of DDA and NMM in 


the engineering fields of rock mass slope instability and rockfall, deformation and failure of underground 


rock mass structure, stress wave propagation and rock blasting, seismic breakage and instability of rock 


)mass structure, and rock mass anchoring. On this basis, the main problems and challenges of the research 


< \and application of DDA and NMM are abstracted and discussed, including the calculation control parame- 


ters, calculation efficiency, multi-physics fields, three-dimensional methods, and the commercial devel- 


*opment of the programs, and corresponding countermeasures are proposed. 


Key words: rock mechanics; discontinuous deformation analysis ( DDA ) ; numerical manifold method 


( NMM) 


岩石 作为 一 种 复杂 地 质 材 料 由 连续 岩石 介质 及 
不 连续 面 如 节理 .断层 和 层 理 等 组 成 ,具有 不 连续 、 
非 均 质 和 各 向 异性 等 特点 。 岩 石 力学 界 发 展 了 各 种 
基于 连续 介质 力学 和 非 连续 介质 力学 的 数值 方法 ， 
用 以 对 岩石 的 各 类 变形 破坏 力学 和 工程 问题 进行 
值 模 拟 。 

连续 数值 方法 如 有 限 单元 法 (finite element 
method, FEM) 有限 差分 法 ( finite difference method , 
FDM) 和 边界 单元 法 (boundary element method, 
BEM) 等 ,在 处 理 复杂 节理 岩 体 以 及 岩石 的 复杂 开裂 
破坏 问题 时 仍 有 很 大 的 局 限 , 且 不 能 很 好 地 处 理 块 
体 旋 转 与 脱离 等 大 位 移 问题 。 为 应 对 上 述 问题 ,可 
进行 连续 - 非 连续 数值 方法 的 耦合 ,如 有 限 元 单元 法 
与 离散 单元 法 (discrete element method, DEM ) Ay #8 


合 … ,或 对 连续 数值 方法 进行 扩展 ,如 扩展 有 限 元 方 
法 (extended finite element method , XFEM ) ?' ,能 较 好 
地 处 理 岩 石 开裂 破坏 问题 。 

在 非 连续 数值 方法 中 ， 离 散 单元 法 ”与 非 连续 
变形 分 析 (discontinuous 
DDA) ”方法 的 基本 思路 都 是 将 分 析 对 象 视 作 离散 
体 的 集合 ， 需 定义 离散 体 间 的 接触 并 在 分 析 过 程 中 
不 断 更 新 接触 信息 ， 适 合 于 模拟 离散 颗粒 或 块 体 的 
力学 行为 。 相 较 于 基于 牛顿 第 二 定律 并 采用 显 式 求 
解 格 式 的 DEM，DDA 则 是 基于 最 小 势能 原理 的 隐 
式 数 值 方 法 ， 具 有 完备 的 运动 学 理论 和 严格 的 平衡 
假定 ， 能 够 在 接触 状态 开 闭 迭代 收敛 的 基础 上 模拟 
可 变形 块 体系 统 的 项 力 和 动力 学 行为 。 

岩石 力学 与 工程 问题 常常 是 连续 与 非 连续 问题 


deformation analysis, 


的 耦合 ,涉及 连续 岩石 介 质 的 变形 破裂 、 不 连续 面 的 
变形 以 及 岩 体 的 大 位 移 等 ,数值 流 形 方法 (numerical 
manifold method,NMM) ”能 很 好 地 适应 此 类 问题 的 
求解 。NMM 采用 数学 和 物理 两 套 履 盖 系 统 , 将 自由 
度 定义 在 物理 覆盖 上 ,实现 了 连续 和 非 连 续 问题 在 
同一 框架 下 的 理论 统一 描述 。NMM 融入 了 FEM 和 
解析 法 的 思想 ,并 同 DDA 一 样 在 接触 开 闭 迄 代 收敛 
的 基础 上 实现 静 力 和 动力 学 问题 的 统一 格式 描述 和 
隐 式 求解 。 当 在 NMM 中 采用 有 限 元 的 网 格 剖 分 来 
生成 数学 覆盖 时 ,NMM 在 形式 上 可 退化 到 FEM; "5 
在 NMM 中 定义 独立 覆盖 描述 离散 体 时 ,NMM 可 退 
化 到 DDA。 因 此 ,FEM 和 DDA 可 以 认为 是 NMM 分 
别 应 用 于 连续 介质 和 离散 体系 统 时 的 特例 。 

XL DDA 和 NMM 分 别 自 上 世纪 80 和 90 年 代 初 提 
HDK RIBERA Jet Re HRN TB 
在 为 学 与 工程 问题 的 重要 方法 。 目 前 , 相 较 于 二 维 
DDA 和 NMM ,由 于 涉及 到 更 为 复杂 的 接触 检测 和 处 
理 等 问题 ,三 维 DDA 和 NMM 的 研究 与 应 用 还 明显 
不 是 。 本 研究 对 DDA 和 NMM 的 理论 与 算法 研究 现 
状态 及 其 在 岩石 力学 与 工程 领域 的 研究 与 应 用 现 
状 洋行 回顾 与 评述 ,在 此 基础 上 总 结 这 2 种 方法 在 
三 帘 与 应 用 上 面临 的 主要 间 题 与 挑战 ,并 提出 应 对 
策略 及 对 未 来 工作 的 展望。 


1 5 理论 与 算法 研究 


-针对 DDA 模拟 面临 的 控制 参数 敏感 . 块 体 旋转 
虚假 体积 膨胀 及 能 量 损 耗 等 问题 ,诸多 学 者 提出 了 
相应 的 优化 方法 。 子 块 体 方法 、 基 于 节点 的 DDA 
( NDDA) 方 法 及 颗粒 DDA 方法 被 提出 以 增强 DDA 
模拟 岩石 变形 和 开裂 破坏 的 能 力 。 高 阶 DDA 及 显 
式 与 并 行 算法 的 研究 提高 了 计算 精度 及 效率 。 在 
NMM 的 研究 中 ,诸多 开裂 算法 的 发 展 提升 了 NMM 
模拟 裂纹 尖端 应 力 奇 异性 、 复 杂 裂 纹 扩展 等 问题 的 
能 力 。 针 对 NMM 独特 的 双重 覆盖 系统 ,有 限 履 盖 、 
独立 覆盖 .部 分 重 等 覆盖 及 覆盖 的 细 化 方法 被 提出 
并 得 到 深入 研究 。 围 绕 单纯 形 积分 及 高 阶 NMM 的 
研究 则 显著 提升 了 计算 精度 。 无 网 格 NMM 及 粒子 
NMM 的 提出 进一步 拓展 了 NMM 的 适用 范围 。 此 
外 ,DDA 和 NMM 中 各 类 材料 本 构 模型 的 租 入 与 发 
展 则 增强 了 这 两 种 方法 在 材料 变形 模拟 中 的 适用 
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性 。 针 对 DDA 和 NMM 在 不 连续 面 显 式 模拟 中 接触 
检测 和 处 理 的 难点 问题 ,起 初 诸多 研究 围绕 接触 检 
测 直 接 法 和 罚 函数 法 开展 ,在 接触 理论 提出 后 ,人 研究 
则 转向 接触 理论 算法 的 实现 和 完善 。 


1.1 DDA 


时 间 步 长 , 步 位 移 比 和 接触 罚 弹 签 刚度 是 DDA 
模拟 的 几 个 关键 计算 控制 参数 。 邬 爱 清 等 加 的 研究 
表明 ,合适 的 时 间 步 长 和 弹簧 刚度 的 取 值 组 合 构成 
一 个 单 连通 参数 取 值 域 。 针 对 DDA 中 块 体 旋 转 时 
的 虚假 体积 膨胀 问题 , Wu SEU 最初 提出 了 接触 后 
调整 法 。 巩 师 林 等 中 则 引入 先 变形 、 后 转动 的 块 体 
位 移 增 量 表达 式 ,考虑 了 块 体 转动 时 离心 力 和 科 氏 
力 引 起 的 变形 。 江 痢 等 ' ”建立 基于 虚 单 元 法 的 
DDA 格式 ,从 块 体 独立 位 移 空 间 上 双 线 性 格式 的 计 
算出 发 ,有 效 克制 了 旋转 膨胀 。 针 对 DDA 块 体 系统 
的 能 量 损耗 问题 ,Jiang 等 “通过 引入 猪 性 阻尼 分 量 
吸收 离散 块 体 的 动能 ,建立 了 考虑 阻尼 效应 的 离散 
块 体系 统 整体 运动 方程 。Yu 等 '“ 则 将 线性 阻尼 器 
施加 在 接触 位 置 ,并 给 出 了 由 给 定 阻尼 比 计算 阻尼 
系数 的 半 解 析 解 。 

一 阶 DDA 的 单个 块 体 具 有 常 应 力 ,针对 岩石 的 
开裂 破坏 模拟 问题 ,Lin 等 "5 提出 了 子 块 体 方法 , 即 
将 连续 岩石 区 域 离散 成 为 由 虚拟 节理 灰 接 的 DDA 
子 块 体 单元 ,从 而 实现 连续 岩石 的 变形 和 破裂 模 拟 。 
Ning 等 (和 Jiao 等 "基于 子 块 体 间 的 接触 应 力 判 
断 沿 子 块 体 单元 边界 的 拉 伸 和 剪 切 开裂 破坏 。 为 减 
小 子 块 体 单元 划分 对 开裂 破坏 模拟 的 影响 ,Ni ES) 
根据 临近 子 块 的 应 力 状 态 判断 沿 虚 拟 节理 面 的 开裂 
破坏 。 张 开 雨 等 "等 则 在 子 块 体 DDA 方法 中 引入 
T Voronoi 多 边 形 离散 形式 。Xia GEO FEF Le pk 
DDA 方法 中 提出 一 种 统一 拉 伸 断裂 模型 用 以 模拟 
各 种 脆性 材料 的 断裂 行为 。 此 外 ,Bao 等 ”提出 了 
NDDA ,该 方法 在 块 体内 部 定义 有 限 元 网 格 ,并 将 网 
格 线 视 为 潜在 裂纹 。 随 后 Tian 等 ”在 NDDA 方法 
中 引入 了 双重 最 小 化 过 程 ,从 而 获得 了 更 准确 的 开 
袭 破坏 模拟 结果 。 在 材料 本 构 模 型 的 嵌 人 与 发 展 方 
ja PEST E DDA 中 将 非 线 性 问题 转化 为 线性 
问题 模拟 块 体 的 物理 非 线性 变形 ,并 建立 了 弹 塑 性 
本 构 关系 。 江 成 “在 DDA 中 实现 了 弹 脆 性 本 构 模 
型 损伤 本 构 模 型 和 Johnson-Cook 本 构 模 型 。 为 描 
述 岩 石 的 非 均 质 特 点 , 杨 正 '” 和 Ning 等 .还 将 基 
于 Weibull 函数 的 岩石 细 观 非 均 质 统计 本 构 模 型 引 
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A DDA ,用 以 模拟 非 均 质 岩石 的 变形 破坏 (图 1) 。 

DDA 本 身 主 要 用 于 模拟 具有 任意 形状 的 多 边 
形 或 多 面体 块 体 , 学 者 们 也 试图 发 展 圆 形 或 球形 颗 
粒 的 DDA 方法 。Ke 460 提出 了 刚性 圆 形 颗粒 
DDA 模型 ,在 此 基础 上 ,Thomas 等 .进一步 提出 了 
圆 形 颗粒 团 簇 模型 以 更 为 真实 的 模拟 粒状 材料 。 郭 
Jt MURS BUD BUR. S TEE E VUE He ood 
合 的 DDA Wik. Jiao EO" 发 展 了 三 维 球 颗粒 DDA 
方法 ,并 建立 了 相 适应 的 黏 结 开裂 算法 ,图 2 为 三 点 
弯曲 实验 的 模拟 结果 。 在 高 阶 DDA 及 显 式 与 并 行 
计算 研究 方面 , Hsiung 等 ”| 在 DDA 中 引入 高 阶 多 
项 式 位 移 函 数 ,使 块 体 复杂 应 力 场 和 应 变 场 的 精确 
模拟 成 为 可 能 。Yang 40^ 在 显 式 时 间 积 分 方案 的 
基础 上 ,提出 了 显 式 DDA 方法 。 张 洪 等 ' 引 详细 讨 
论 子 基于 中 心 差分 法 的 三 维 显 式 DDA 格式 。Peng 
竺 将 并 行 块 雅 可 比方 法 和 预 处 理 共 亏 梯度 侈 代 
求解 方法 实现 到 三 维 DDA 中 ,随后 进一步 提出 了 基 
OpenMP 的 三 维 DDA 全 阶段 并 行 化 方法 055 。 王 
等 5 实现 了 三 维 DDA 总 体 平衡 方程 组 的 雅 可 
比 巴 处 理 共 罗 梯度 法 的 CUDA 并 行 求解 。 


图 1 预制 裂纹 非 均 质 岩石 圆 盘 径 
向 压缩 破坏 的 DDA 模拟 


TYig-1 DDA simulation of heterogeneous rock disk fracturing 
[25] 


x 
G 
= 

E 


Oo with pre-exiting crack under radial compression 


图 2 三 点 弯曲 破坏 的 三 维 球 颗粒 DDA 模拟 ” 
Fig.2 3D spherical particle DDA simulation 
[30] 


of three-point bending failure 


1.2 NMM 


NMM 连续 - 非 连续 相 统 一 的 特点 让 其 在 岩石 的 
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开裂 破坏 模拟 上 具有 显著 优势 。Ma A807 和 Zhang 
等 中 针对 裂纹 尖端 的 奇异 性 引入 奇异 物理 覆盖 , 利 
用 NMM 模拟 了 复杂 多 支 及 交叉 裂纹 问题 。 徐 栋 栋 
等 (在 物理 覆 羡 上 采用 位 移 函 数 的 一 阶 泰勒 展开 ， 
并 在 裂纹 尖端 附近 的 物理 覆盖 上 增加 扩充 位 移 函 
数 ,以 更 好 模拟 裂纹 尖端 的 应 力 奇异 性 。Chiou 
等 后 提出 了 结合 虚拟 裂纹 扩展 方法 的 NMM, 采 用 
最 大 切 应 力 准 则 确定 裂纹 扩展 方向 ,模拟 了 混合 型 
Bory Ew, Ning 等 “在 NMM 中 提出 了 一 种 基 
于 拉 伸 强度 准则 和 Mohr-Coulomb 准则 的 开裂 算法 ， 
并 模拟 了 下 盘 边 坡 的 破裂 失 稳 ;在 此 基础 上 , Kang 
等 (1 对 和 岩石 试 样 的 单 裂纹 、 双 裂纹 和 多 裂纹 扩展 和 
贯通 问题 进行 了 系统 性 模拟 ,展示 了 结果 的 正确 性 
和 NMM 对 复杂 开裂 破坏 问题 的 处 理 能 力 (图 3)。 
针对 动态 裂纹 扩展 ,Zheng 等 “将 时 间 离 散 置 于 空 
间 离 散 之 前 ,解决 了 裂纹 扩展 中 自由 度 在 单 时 步 内 
不 断 变 化 和 能 量 守恒 的 问题 。 针 对 三 维 裂纹 问题 ， 
Yang 45^" 基于 最 大 拉 应 力 准 则 开发 了 模拟 三 维 裂 
颖 扩展 的 NMM, 


图 3 预制 裂纹 岩石 试 样 单 轴 压 缩 开 裂 
破坏 的 NMM Epp 7 


Fig.3 NMM simulation of rock specimen fracturing 
[42] 


with pre-existing cracks under uniaxial compression 

双重 覆盖 系统 是 NMM 实现 连续 - 非 连续 相 统一 
的 核心 和 基础 。 二 维 NMM 多 通过 规则 三 角形 网 格 
或 有 限 元 三 角形 网 格 来 生成 有 限 覆 盖 系 统 。 在 窗 盖 
系统 研究 方面 , Yang 等 “ 研究 了 数学 覆盖 的 细 化 问 
题 ,详细 介绍 了 数学 覆盖 在 不 同 边界 条 件 下 的 更 新 
过 程 。Cai 等 “提出 了 一 种 简洁 的 广义 覆盖 生成 方 
法 ,适用 于 复杂 不 连续 问题 的 处 理 。 祁 勇 峰 等 ”| 提 
出 了 部 分 重 欠 履 盖 的 NMM ,解决 了 禾 盖 不 匹配 导致 
的 结构 关键 部 位 计算 精度 下 降 的 问题 。 武 奢 等 ” 研 
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究 了 针对 复杂 岩 体 力学 问题 的 有 限 上 覆盖 与 独立 覆盖 
混合 使 用 方法 ,模拟 分 析 了 复杂 结构 岩 质 边 坡 问题 
(图 4) 。 此 外 ,NMM 通过 单纯 形 积分 实现 任意 形状 
流 形 单元 上 的 精确 积分 运算 。 在 单纯 形 积 分 的 研究 
方面 ,Kourepinis 等 “提出 了 一 种 对 任意 阶 覆 盖 位 
移 函 数 进行 单纯 形 积分 的 方法 , Liu 等 ”推导 了 二 
阶 位 移 函 数 的 单纯 形 积分 表达 式 。 


BLAIR 


图 4 复杂 岩 质 边 坡 失 稳 破坏 的 NMM 模拟 
ENI Fig.4 NMM simulation of instability and breakage 
d [8] 
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of a complex rock slope 
-一 在 高 阶 NMM 的 研究 方面 ,An ^60" 提出 了 基于 
分 解 的 高 阶 数值 方法 的 系统 方程 秩 亏 问题 , 建 
立 于 网 格 秩 亏 数 的 定量 预测 方法 。Xu 等 5 THE 
Wis NMM 中 的 线性 相关 问题 ,基于 FEM 三 角形 
板 欠 元 推导 出 一 种 新 的 局 部 逼近 。 苏 海 东 等 提 
出 了 应 用 软件 自动 推导 公式 和 生成 程序 代码 的 简便 
方法 ,开发 了 高 阶 NMM 的 二 维和 三 维 静 力 分 析 程 
序 。 苏 海 东 等 3 提出 了 高 阶 NMM 的 两 种 初 应 力 处 
理 方 法 。Yang 等 .通过 四 边 形 网 格 形成 数学 覆 
盖 , 提 出 了 一 种 基于 约束 和 正 交 最 小 二 乘 方法 局 部 
逼近 的 高 阶 NMM。 在 材料 本 构 模 型 的 嵌入 与 发 展 
方面 ,Zeng A 05! 将 增 量 内 时 理论 推广 到 NMM 中 ， 
提出 了 一 种 弹 塑 性 分 析 的 内 时 化 算法 。He 等 在 
NMM 中 引入 基于 广义 Kelvin-Voigt EA HA Fk HEIM 
量 本 构 关系 ,以 更 好 地 模拟 黏 弹性 材料 的 蠕 变 裂 纹 
WMA. Kang 等 ”在 NMM 中 引入 黏 弹 本 构 模 型 、 
弹 秋 塑 性 本 构 模 型 描述 高 聚 物 黏 结 炸药 中 聚合 物 秋 
结 剂 ,炸药 颗粒 的 变形 。 

在 其 他 形式 的 NMM 研究 方面 ,Li 等 基于 单 
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[LA A AAT RA UE HEL FICS NMM, 
并 随后 研究 了 无 网 格 NMM 中 裂纹 尖端 数值 精度 的 
提高 方法 '” 。Gao 等 ' 基于 复 变 量 移动 最 小 二 乘 
逼近 和 有 限 履 盖 理 论 ,提出 复 变量 无 网 格 NMM , 并 
从 最 小 势能 原理 出 发 给 出 了 适用 于 断裂 问题 的 相应 
公式 。 李 伟 等 ' “通过 整合 裂纹 尖端 附近 的 数学 覆 
盖 生 成 复合 覆盖 , 并 在 其 上 定义 有 限 项 的 位 移 Wil- 
liams 级 数 ,建立 了 求解 线 弹 性 断裂 问题 的 无 网 格 
NMM。 此 外 ,Sun 等 !” 还 通过 数学 覆盖 系统 来 描述 
基于 粒子 的 运动 和 变形 场 , 提 出 了 粒子 NMM, Li 
等 (中 则 进一步 发 展 该 方法 用 于 岩石 动态 破裂 行为 
的 模拟 (图 5)。Wu 等 ' 引 通过 岩石 的 Voronoi 微观 
颗粒 划分 实现 了 一 种 岩石 微观 力学 特性 NMM 模拟 
技术 。Zhou 等 ' 外 则 提出 一 种 基于 Voronoi 颗粒 的 改 
进 NMM 来 模拟 细 观 尺度 下 完整 岩石 的 变形 破坏 。 
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图 5 炮 孔 爆破 的 粒子 NMM 模拟 ~ 
Fig.5 Particle NMM simulation of borehole blasting 1 


1.3 # f 


不 连续 面 及 离散 体 间 的 接触 是 非 连续 计算 的 重 
要 特征 。 接 触 的 检测 及 处 理 是 DDA 和 NMM 发 展 过 
程 中 面临 的 最 困难 的 一 类 问题 “1! 。 围 绕 二 维 接触 ， 
Keneti 等 ' 引 提出 一 种 凸 块 接触 检测 方法 ,通过 将 两 
个 逼近 面 中 的 一 个 视 为 “ 主 平面 "并 搜索 一 定 范围 
内 另 一 个 面 的 顶点 ,在 此 基础 上 识别 相 邻 块 间 的 主 
要 接触 位 置 及 类 型 。Zhuang 45/9" 提出 一 种 基于 履 
盖 的 不 规则 凸 多 边 形 接触 检测 策略 ,利用 顶点 - 边 和 
边 -顶点 接触 覆盖 描述 多 边 形 间 的 接触 约束 。 葛 发 
Ik 7" 提出 了 刚性 圆 盘 与 边 . 凸 角 和 四 角 的 接触 检 
测 方法 。 对 于 三 维 接触 , Wu 等 "在 三 维 DDA 中 利 
用 多 重 覆 盖 系 统 和 几何 子 单元 分 解 将 邻 块 的 搜索 与 
接触 类 型 的 判定 联系 起 来 ,显著 减少 了 接触 检测 量 
和 迭代 次 数 。Zhang 等 "针对 三 维 任意 形状 多 面 
体 ,通过 识别 多 面体 几何 元 素 的 一 般 特征 和 关系 , 利 
用 邻 域 搜索 、 接 触 类 型 检验 和 侵入 识别 进行 接触 的 
检测 。He 7) WEE RE = Ae NMM 中 的 接触 检测 问 
题 提 出 一 种 分 层 接触 系统 ,实现 了 三 维 非 网 格 模型 
的 一 般 接触 求解 过 程 。 
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对 于 接触 的 处 理 ,DDA 和 NMM 都 主要 采用 罚 
函数 法 。Jiao 等 "在 DDA 中 提出 了 一 种 两 阶段 接 
触 模型 来 模拟 节理 岩 体 ，Gong 等 中 也 将 两 阶段 模 
型 引入 DDA 中 描述 边 - 边 接触 行为 ,以 模拟 岩 体 的 
拉 伸 破坏 过 程 。Zheng 等 1e DDA 子 块 体 单元 的 
分 布 茜 结 模拟 方法 中 ,考虑 了 线性 及 非 线 性 条 结 模 
型 。Liu 等 "采用 BB 和 双 曲 节理 模型 分 别 对 罚 弹 
得 的 法 向 和 前 切 刚度 进行 修正 ,在 NMM 中 提出 了 
一 种 描述 节理 力学 行为 的 修正 接触 模型 。Kang 
等 3 在 NMM 中 引入 双 线性 狐 接 接触 模型 描述 高 
聚 物 黏 结 炸药 中 颗粒 与 黏 结 剂 界面 的 变形 。 对 于 接 
触 的 开 闭 迭代 ,Jiang 等 ”在 导出 三 维 DDA 点 - 面 接 
触 的 罚 弹 簧 和 摩擦 力 子 和 矩阵 的 基础 上 ,给 出 了 不 同 
接触 状态 的 开 闭 迭代 准则 。Zhou 等 ” 则 在 NMM 
中 通过 删除 接触 边界 上 的 中 间 点 来 实现 更 高 的 接触 
mou a 
的 通 取 非常 敏感 ,一 些 研究 考虑 了 对 罚 函 数 法 的 改 
进 5 其 至 考虑 取缔 罚 弹 簧 。Amadei 45" E DDA 中 
采用 增 广 拉 格 朗 日 乘 子 法 代替 罚 函 数 法 ,弱化 了 罚 
SERIU. Zheng 等 ”以 接触 力作 为 基本 变 
量 Y 应 用 有 限 维 变 分 不 等 式 理论 建立 DDA 的 对 侦 形 
式 ,达到 了 完全 移 除 虚拟 罚 弹 签 的 目的 。 

OJEE DDA 和 NMM 以 及 后 续 众 多 学 者 的 研究 
中 5 接触 检测 都 是 采用 直接 法 ,尤其 是 针对 三 维 问 

是 可 能 面临 接触 检测 的 失败 。 关 于 接触 检测 
fiae cV EE Té Shi'” ”提出 了 适用 于 任意 二 维和 
三 重 几 何 体 间接 触 检测 的 接触 理论 ,也 称 为 进入 块 
体 理 论 。 该 理论 通过 定义 接触 覆盖 并 形成 进入 块 ， 
进 希 将 任意 两 个 块 体 4 和 8 间 的 接触 关系 简化 为 参 
考点 a6 与 进入 块 E(4, B) 间 的 关系 ,其 中 (4,B) = 
B-A+ao。 图 6(a) 为 2 个 二 维 多 边 形 块 体 的 E(4， 
B) 示 意图 。 接 触 理论 极 大 地 精简 了 接触 检测 的 计 
算 ,更 从 理论 上 实现 了 任意 形状 块 体 间接 触 检 测 的 
完备 性 ,从 而 突破 了 三 维 DDA 和 NMM 发 展 面临 的 
关键 障碍 。 在 此 基础 上 , Lin 等 提出 了 在 三 维 
DDA 程序 中 引入 接触 理论 算法 的 粗略 策略 。Ni 
等 (中 1 则 详细 提出 了 二 维 任意 多 边 形 块 体 和 三 维 凸 
多 面体 块 体 的 接触 理论 算法 并 在 DDA 中 进行 了 实 
现 。 在 该 算法 中 ,多 面体 块 体 间 的 六 种 基本 接触 类 
型 被 处 理 为 点 - 面 和 非 平行 边 - 边 两 种 接触 类 型 ,并 
同时 对 点 -点 和 点 - 边 接触 进行 特殊 处 理 ,保证 了 算 
法 的 有 效 性 和 效率 。 图 6(b ) 为 基于 接触 理论 的 三 
维 块 体 群 滑动 的 DDA 模拟 。 
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(b) 三 维 块 体 群 滑动 模拟 下 
图 6 接触 理论 及 其 在 三 维 DDA 中 的 应 用 


Fig.6 Contact theory and its application in 3D DDA 


2 ”岩石 力学 与 工程 应 用 研究 


岩 体 边 坡 失 稳 与 落石 是 常见 的 地 质 灾害 ,二 维 
DDA 常 被 用 来 模拟 分 析 边 坡 的 稳定 性 及 失 稳 破 坏 ， 
针对 落石 冲击 多 通过 三 维 DDA 进行 模拟 ,NMM 则 
常用 来 模拟 边 坡 的 渐进 破坏 。 在 隧道 及 地 下 空间 工 
程 中 ,DDA 与 NMM 被 用 来 模拟 分 析 隧 道 及 砚 室 围 
宕 的 应 力 分 布 与 变形 规律 以 及 破坏 失 稳 。 为 更 好 地 
模拟 岩石 动力 学 问题 中 涉及 的 应 力 波 传播 问题 ,于 
绕 DDA 与 NMM 的 动力 学 边界 条 件 模型 开展 了 大 量 
研究 ,在 此 基础 上 在 DDA 中 引入 爆炸 应 力 波 与 爆 生 
气体 的 作用 模拟 分 析 岩 石 爆破 问题 。 针 对 岩 体 结构 
的 地 震 响应 这 一 类 重要 岩石 动力 学 问题 , DDA 及 
NMM 中 的 地 震 加 载 方式 得 到 了 较 多 人 研究, 地震 诱发 
的 岩 体 滑坡 及 地 下 岩 体 结构 失 稳 灾害 得 到 了 较 好 的 
模拟 分 析 。 和 针对 锚 杆 加 固 这 一 重要 岩 体 支 护 手段 ， 
DDA 与 NMM 主要 围绕 二 维 锚固 模型 的 开发 进行 研 
FE ,并 通过 拉 拔 与 锚固 实验 进行 验证 。 


2.1 岩 体 边 坡 失 稳 与 落石 


在 岩 体 边 坡 失 稳 的 DDA 模拟 方面 ,Sitar 459 
通过 对 两 个 典型 边 坡 破 坏 实 例 的 模拟 ,分 析 了 结构 


面 位 置 和 数量 对 边 坡 稳定 性 预测 的 显著 影响 。Chen 
A5 5 模拟 了 新 磨 边 坡 的 失 稳 后 行为 ,DDA 确定 的 破 
坏 过 程 与 通过 地 震 信号 分 析 得 到 的 结果 有 较 好 的 一 
致 性 。 针 对 边 坡 倾倒 失 稳 , 孙 东亚 等 中 用 DDA 对 
基于 极限 平衡 原理 的 分 析 结果 进行 验证 ,研究 了 全 
倒 破 坏 的 失 稳 变形 机 理 。 涉 及 水 利 工程 时 , 沈 振 中 
等 考虑 水 库 水 位 又 降 时 地 下 水 渗透 压力 的 作用 ， 
通过 模拟 得 到 了 多 种 水 位 变化 条 件 下 岩 体 边 坡 的 稳 
[1233 2851 7 DIE 在 自动 强度 折 减 法 的 基 
础 上 ,分 析 了 水 库 岩 体 边 坡 的 稳定 性 并 得 到 岩 体 边 
坡 安全 系数 随 水 位 的 变化 规律 。 

在 落石 的 DDA 模拟 方面 , Chen 257 开展 了 滚 
石 过 程 的 三 维 DDA 模拟 (图 7) ,并 对 比 了 二 维和 三 
维 DDA 预测 落石 轨迹 和 动力 学 行为 的 效果 。Ma 
等 灸 ;开发 了 双 目 立体 视觉 野外 实验 系统 并 结合 三 
ADDA 模拟 结果 ,定量 得 到 了 不 同 条 件 下 落石 运动 
m 1i 
ETRAS" 以 西藏 自治 区 G318 国道 KA580 典型 滑坡 
将 着 景 的 研究 表明 , 三维 DDA 能 够 有 效 模拟 巨石 项 
ee 
| 


力学 过 程 。 
= ® Case | 
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图 7 落石 深 动 三 维 DDA a”! 
Fig.7 3D DDA simulation of rockfall rolling: 
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在 岩 体 边 坡 失 稳 的 NMM 模拟 方面 , Yang 5 
利用 基于 多 层 数学 覆盖 的 改进 NMM 对 土石 混合 体 
边 坡 进 行 了 稳定 性 分 析 。Wong 等 .利用 NMM 成 
功 模 拟 了 兰 体 边 坡 渐进 破坏 中 的 袭 纹 萌 生 、 扩 展 、 合 
并 、 退 化 直至 最 终 的 灾难 性 破坏 全 过 程 。He SEO) 
提出 了 分 析 节 理 岩 体 边 坡 稳定 性 的 三 维 NMM ,讨论 
了 其 在 岩 质 边 坡 稳定 性 分 析 中 存在 的 主要 问题 , 强 
调 了 NMM 在 边 坡 稳定 性 分 析 中 的 突出 特点 。Liu 
等 5%] 利 用 三 维 独 立 覆 盖 NMM(ICMM3D ) 计算 边 坡 
的 应 力 场 ,利用 矢量 和 法 (VSM ) 计算 安全 系数 ,有 效 
开展 了 三 维 边 坡 的 稳定 性 计算 。 针 对 边 坡 的 倾倒 破 
坏 , 张 国 新 等 “结合 极限 平衡 原理 和 NMM 计算 ,得 
出 了 单个 块 体 及 整体 结构 的 安全 系数 。 王 欢 欢 
等 "以 青藏 高 原 金沙 江 流域 的 西藏 昌都 地 区 贡 扎 
倾倒 滑坡 为 例 , 利 用 NMM 研究 了 重力 作用 下 反 倾 
层 状 岩 体 边 坡 在 不 同 岩 层 倾角 下 的 破坏 特征 。 此 
外 ,Wu 等 ' 呈 还 在 NMM 中 结合 流动 模型 和 裂隙 演化 
技术 ,模拟 了 岩 体 边 坡 在 暴雨 作用 下 的 破坏 过 程 。 


2.2 地 下 岩 体 结构 的 变形 与 破坏 


地 下 开 挖 工程 会 破坏 岩 体 原 有 的 平衡 状态 ,从 
而 导致 围 岩 体 结构 的 变形 甚至 破坏 和 失 稳 。 采 用 
DDA 方法 ,Wu 等 ”模拟 分 析 了 隧道 开 挖 过 程 中 崇 
体 的 应 力 拱 效 应 ,揭示 了 室内 实验 中 垂直 应 力 和 地 
表 沉降 不 对 称 的 原因 。 刘 君 等 "模拟 了 节理 岩 体 
中 的 隧道 开 控 过 程 ,研究 了 不 同 节理 倾角 条 件 下 兰 
体 的 应 力 分 布 特性 以 及 开 挖 后 隧道 围 岩 的 变形 和 应 
力 分 布 规律 。Bakun-Mazor 等 "对 以 色 列 耶路撒冷 
旧 城 下 层 状 节理 岩 体 中 的 一 个 独立 无 支 护 地 下 采 石 
场 进 行 模拟 ( 图 8) ,结果 表明 随 着 节理 长 度 的 增加 
和 涯 桥 长 度 的 减 小 , 岩 体 的 竖 向 变形 失 稳 增 大 。 
Yeung 等 "| 研究 了 不 同 埋 深 \、 节 理 方位 节理 间距 
和 节理 摩擦 角 组 合 情 况 对 埠 道 稳定 性 的 影响 。Kim 
等 ("进行 了 韩国 某 实 际 隧道 的 开 挖 模拟 ,研究 了 询 
隙 流动 、 开 挖 顺序 和 加 固 措 施 对 隧道 稳定 性 的 影响 。 
Tsesarsky 等 ("通过 模拟 重点 研究 了 节理 间距 和 摩 
擦 角 对 地 下 砚 室 开 挖 稳定 性 的 影响 。 乌 爱 清 等 
以 清江 水 布 十 水 电工 程 地 下 厂房 为 例 , 用 DDA 对 复 
林地 质 条 件 下 地 下 厂房 围 岩 的 变形 与 破坏 特征 开展 
了 研究 ,重点 分 析 了 地 应 力 水 平 、 锚 固 、 岩 体 结构 及 
结构 面 强度 的 影响 。 采 用 NMM 方法 ,Tal 等 ”成 
功 模 拟 了 高 地 应 力 条 件 下 非 连 续 岩 体 地 下 硼 室 的 稳 
定性 问题 。 曹 文 贵 等 ”在 NMM 中 提出 了 公路 路 
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基 兰 溶 项 板 稳定 性 的 评价 方法 ,分 析 了 岩溶 顶板 稳 
定性 与 岩溶 埋 深 和 溶洞 几何 形态 等 因素 的 关系 。 


L-5, B=1 


L-10, B=2 


图 8 地 下 采 石 场 上 部 围 岩 失 稳 的 DDA 模拟 
Fig.8 DDA simulation of roof rock instability 
[101] 
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of an underground quarry 
e 
20 应 力 波 传播 与 岩石 爆破 


pene! 波 在 人 工 截 取 边 界 处 的 反射 问题 ,Bao 
等 WW 发展 了 一 种 对 纵波 和 横 波 都 有 较 好 吸收 能 

EDDA 黏 性 边界 条 件 。Wa 等 ”在 NMM 中 提出 
ON 种 改进 的 一 阶 Higdon 吸收 边界 。 在 应 力 波 传 
播 的 模拟 分 析 方 面 ,Gu 等 ("通过 DDA 研究 了 块 体 

.边界 条 件 和 入 射 波 频率 对 波 传播 的 影响 ,分 析 
了 义 响 波 传播 模拟 精度 的 因素 。Ning 等 "建立 了 
DD 接触 罚 弹 得 刚度 与 节理 刚度 的 对 应 关系 模型 ， 
通 时 模拟 分 析 了 节理 刚度 、 节 理 强 度 以 及 载荷 频率 
与 见 射 角 等 因素 对 节理 岩 体 应 力 波 传播 的 影响 。 
Fan 等 5221 采用 NMM 模拟 分 析 了 应 力 波 穿越 裂 姻 时 
的 频率 依赖 关系 ,讨论 了 裂缝 刚度 对 波 衰减 和 有 效 
波 速 的 影响 。Zhou 460 进而 研究 了 波 在 微 裂隙 岩 
体 中 的 传播 规律 ,NMM 模拟 表明 微 裂缝 长 度 、 数 量 
等 几何 参数 对 有 效 速度 .衰减 速率 .衰减 系数 和 波 数 
有 显著 影响 。Zhao an] 通过 引入 Newmark 方程 组 
解决 NMM 原 系统 方程 导致 的 人 工 数值 阻尼 ,并 对 
接触 方案 进行 优化 以 更 好 模拟 应 力 波 通过 不 连续 面 
的 传播 问题 。Fan 等 25 采用 NMM 研究 了 裂纹 对 应 
力 波 持续 时 间 和 振幅 的 显著 影响 特性 。 

在 岩石 爆破 的 模拟 方面 , Mortazavi 等 "在 
DDA 中 发 展 了 炮 孔 扩张 模型 ,考虑 了 炮 孔 几 何 参 
数 .完整 岩石 性 质 、 既 有 结构 面 的 几何 分 布 与 力学 性 
质 以 及 炮 孔 压力 对 岩 休 破坏、 抛掷 和 堆积 过 程 的 影 
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响 。Ning 等 “通过 在 DDA 中 引入 子 块 体 单元 
开裂 破坏 模拟 方法 并 考虑 爆 生 气体 在 岩 体 裂缝 中 的 
渗透 加 载 作 用 ,完整 再 现 了 漏斗 和 台阶 抛掷 爆破 中 
的 爆 腔 扩 张 及 岩石 破裂 抛掷 与 堆积 全 过 程 (图 9) 。 
在 后 续 工 作 中 , 郭 双 等 9 还 研究 了 爆炸 应 力 波 和 爆 
生气 体 压 力 共同 作用 下 的 爆破 破 岩 过 程 ,并 考虑 了 
地 应 力 的 影响 。 


图 9 炮 孔 抛掷 爆破 的 DDA 模拟 
Fig.9 DDA simulation of borehole throw blasting ^ 


] 


2.4 岩 体 结构 地 震 破 坏 与 失 稳 


Wu 45 介绍 了 DDA 中 地 震 载荷 的 其 中 3 种 
施加 方式 ,分 别 是 在 滑 体 上 加 载 加 速度 时 程 、 在 基 岩 
上 加 载 加 速度 时 程 和 在 基 崖 上 加 载 位 移 时 程 。Ning 
等 (31 详尽 推导 了 地 震 作 用 下 二 维 单 块 体 滑动 的 理 
论 解 ,并 深入 验证 和 分 析 了 单 块 和 多 块 地 震 滑动 的 
DDA 模拟 。 在 地 震 滑 坡 的 模拟 方面 , Hatzor 等 
以 加 速度 时 程 作 为 输入 并 引入 2% 的 动态 阻尼 , 对 
以 色 列 King Herod’s Palace 上 部 阶地 的 高 度 非 连续 
岩 体 边 坡 进行 了 二 维稳 定性 分 析 , 得 到 了 与 实际 情 
况 吻 合 的 模拟 结果 。Wnu 等 "采用 DDA 模拟 研究 
了 地 震 诱 发 的 草 岭 滑坡 中 滑坡 体 的 运动 特性 。 
Zhang 等 "通过 DDA 模拟 研究 了 大 光 包 滑坡 (图 
10) ,揭示 了 近 断 层 地 震 力 对 滑坡 过 程 、 滑 动 距离 和 
破坏 形态 的 显著 影响 。Wang 4U 用 DDA 对 2008 
年 汶川 地 震 诱 发 的 东 河 口 高 速 远程 滑坡 现象 进行 了 
模拟 分 析 , 描述 了 滑坡 的 整个 动力 学 过 程 。Zhang 
等 (在 DDA 中 通过 设置 竺 性 边界 条 件 并 采用 应 力 
输入 方式 施加 地 震 载 荷 , 对 大 岗 山 水 电站 地 下 厂房 
的 地 震 响 应 进行 了 模拟 分 析 。 在 NMM rp, Wei 
Ag 08 开发 了 经 典籍 性 边界 条 件 、 自 由 场 边界 条 件 和 


静 动 态 统一 边界 ,扩展 了 采用 NMM 3EG HAAS 
地 震 分 析 的 能 力 。Yang 等 "” 在 对 NMM 进行 一 系 
列 动 态 边界 条 件 开发 的 基础 上 ,实现 了 岩 体 地 震 响 
应 的 有 效 模拟 。 


x< 图 10 大 光 包 滑坡 的 DDA 模拟 
(Fig.10 DDA simulation of Daguangbao landslide E? 


c 
2:8- ath 8 El 


中 在 DDA 中 ,Nie 5? See Lah op T Acn T 
的 线 弹性 应 变 硬化 模拟 模型 ,并 开展 了 拉 氢 荷载 作 
用 下 锚 杆 系统 的 承载 和 变形 能 力 模拟 分 析 。Ma 
等 ”提出 错 杆 - 岩 体 界面 的 三 线性 稿 结 滑 移 模 型 ， 
通过 锚 杆 拉 拔 实验 研究 了 界面 切 应 力 的 变化 规律 及 
脱 黏 现象 。 寓 尤 军 等 5 采用 多 线性 黏 结 滑 移 界面 
模型 详细 探讨 了 锚 杆 的 滑 移 脱 秋 行 为 (图 11) Nie 
等 20 还 通过 锚 杆 及 岩 体 的 DDA 实体 单元 建 模 , 模 
拟 分 析 了 连续 力学 耦合 ( CMC ) 错 杆 单元 的 黏 接 滑 
移 破坏 行为 。He 等 ("通过 将 错 杆 离散 为 欧 拉 伯 努 
利 染 单 元 ,在 DDA 中 提出 了 一 种 统一 错 杆 模型 。 此 
外 , 张 秀 丽 等 早期 在 DDA 中 建立 了 简单 的 预 应 
力 错 杆 模型 ,模拟 结果 体现 了 预 应 力 错 杆 对 节理 岩 
体 的 加 固 作 用 。 陈 云 娟 等 "在 DDA 中 建立 了 全 长 
前 切 钳 杆 模型 ,并 对 一 公路 隧道 锚固 稳定 性 进行 了 
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分 析 。 姜 清 辉 等 较 早 提出 了 三 维 DDA 中 锚 杆 的 
简化 模型 ,并 通过 算 例 证 明 模 型 能 较 好 地 反映 锚 杆 
的 加 固 效 果 。 在 NMM rp, Wei 等 开发 了 一 种 基 
于 梁 单 元 的 销 杆 模型 ,并 对 水 布 址 水 电站 地 下 厂房 
的 锚 杆 加 固 效果 进行 了 模拟 和 评价 。 董 志 安 等 呈 " 
较 早 也 在 NMM 中 开发 了 锚固 支 护 功 能 ,并 对 某 层 
状 结构 岩 体 地 下 三 室 围 岩 开 挖 变形 问题 进行 了 模拟 
研究 。 
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图 11 锚 杆 拉 拔 DDA 模拟 与 实验 结果 的 比较 
Fig. 11 Comparisons between DDA and experimental results 


of bolt pullout '^ 


3 挑战 与 建议 


DDA 与 NMM Bj fé thi DOK, TE JT RUBUS .接触 检 
W 流 形 柳 盖 三 维 方法 等 方面 开展 了 较 多 的 理论 与 
算法 研究 工作 ,并 在 边 坡 及 地 下 岩 体 结构 失 稳 、 岩 体 
爆破 及 地 震动 力学 响应 、 岩 体 锚 固 等 岩石 力学 与 工 
程 问题 中 得 到 了 广泛 成 功 的 应 用 。 尽 管 取得 了 显著 
成 效 ,但 归纳 起 来 ,笔者 认为 DDA 和 NMM 在 研究 与 
应 用 上 还 面临 如 下 主要 挑战 。 
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1) 控制 参数 问题 

DDA 与 NMM 的 模拟 结果 对 时 间 步 长 . 步 位 移 
比 和 接触 罚 弹簧 刚度 这 3 个 计算 控制 参数 非常 敏 
感 。 尤 其 是 对 于 较为 复杂 的 岩石 力学 问题 ,用 户 需 
建立 对 关键 控制 参数 意义 及 其 相互 关系 的 深刻 理 
解 ,并 积累 丰富 的 使 用 经 验 才能 得 到 可 靠 的 模拟 分 
析 结 果 。DDA 和 NMM 均 采 用 动态 系数 ( 质点 速度 
逐 时 步 继承 比例 ) 来 反映 问题 的 动态 程度 ,其 取 值 一 
般 也 较 难 确定 。 建 议 在 深入 揭示 控制 参数 物理 意义 
的 基础 上 ,建立 更 为 完善 的 参数 自动 确定 方法 ;对 于 
动力 学 计算 ,在 通过 非 线性 材料 和 接触 模型 等 渠道 
进行 与 实际 物理 过 程 相 一 致 的 能 量 消耗 的 基础 上 ， 
可 以 考虑 对 动态 系数 参数 予以 取消 。 此 外 ,接触 罚 
弹 答 在 静 力 和 动力 学 问题 中 具有 截然 不 同 的 意义 ， 
因 邮 也 应 予以 深入 研究 。 
) 计 算 效 率 问题 
DDA 与 NMM 都 是 基于 最 小 势能 原理 的 隐 式 方 
AER 每 个 时 步 都 要 达到 接触 开 闭 迭代 状态 的 收敛 ， 
等 二 个 迭代 步 都 要 进行 系统 矩阵 的 组 装 和 系统 方程 
的 求解 ,再 加 上 每 个 时 步 初始 的 接触 检测 ,因此 在 实 
际 遍 用 过 程 中 ,尤其 是 对 于 大 规模 问题 ,往往 带 来 较 
难 若 受 的 计算 量 。 尽 管 学 者 们 在 显 式 与 并 行 算法 开 
发 -接触 检测 算法 优化 等 方面 开展 了 较 多 的 研究 工 
人 得 隐 式 格式 基础 上 的 并 行 在 效率 提升 的 空间 上 
极为 有 限 ,而 显 式 计算 又 不 利于 保留 接触 开 闭 迭代 
收 令 这 一 重要 优势 和 特征 。 因 此 在 计算 效率 的 提升 
方面 ,建议 在 接触 检测 的 策略 优化 和 并 行 处 理 , 以 及 
接触 开 闭 迭 代 与 显 式 并 行 的 相 容 性 方面 开展 深入 研 
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至 今 仍 非常 不 足 。 尽 管 接触 理论 的 提出 为 扫除 接触 
检测 这 一 重大 发 展 障碍 提供 了 极 具 意义 的 理论 文 
撑 , 但 其 具有 和 鲁 棒 性 的 算法 实现 以 及 在 此 基础 上 的 
接触 处 理 都 还 是 三 维 DDA 和 NMM 发 展 路 上 所 需 解 
决 的 重大 问题 。 由 于 DDA 和 NMM 都 是 对 非 连续 面 
的 显 式 模 拟 ,三 维 模型 的 有 效 和 快速 构建 至 关 重 要 。 
目前 主要 采用 的 块 体 切割 建 模 方 法 多 面临 输入 参数 
量 大 、 算 法 稳定 性 差 等 问题 。 因 此 在 三 维 方法 的 发 
展 方面 ,建议 加 强 基 于 接触 理论 的 接触 检测 算法 研 
究 , 丰 富 三 维 模型 的 建立 手段 ,如 建立 三 维 DDA/ 
NMM 结构 模型 与 主流 三 维 地 质 模型 间 的 联系 等 ,并 
在 计算 效率 的 提升 方面 加 大 研究 力度 。 

5 ) 程 序 的 商业 化 开发 问题 

内 及 国际 关于 DDA 和 NMM 的 研究 基本 上 都 
是 在 石 根 华 先生 的 源 代码 基础 上 各 自 开 发 发 展 , 产 
生 了 较 多 的 重复 性 工作 ,而 更 重要 的 是 无 法 形成 对 
DDA 和 NMM 具有 重要 推广 和 应 用 意义 的 商业 化 程 
序 , 加 之 计算 控制 参数 与 计算 效率 的 制约 ,严重 影响 
了 方法 的 发 展 与 应 用 进程 。 建 议 学 者 及 团队 间 增 强 
交流 ,分 工 合作 ,通过 协同 开展 重大 项 目 申 报 与 实施 
等 ,逐步 形成 成 熟 并 具有 影响 力 的 DDA/NMM 岩石 
力学 商业 软件 ,解决 更 多 具有 重大 意义 的 工程 实际 


问题 。 


4 结论 及 展望 


DDA 与 NMM 以 其 完备 的 理论 实现 了 对 非 连 续 
以 及 连续 - 非 连续 相 耦 合 问题 的 有 效 求解 , 现 已 成 为 


究 $ 对 接触 处 理 方法 的 优化 ,甚至 考虑 对 接触 罚 弹 答 
的 取缔 ,可 能 是 值得 研究 的 重要 方向 。 

3) 多 物理 场 问题 

DDA 与 NMM 中 多 物理 场 耦 合 功能 的 开发 符合 
我 国 次 地 岩石 工程 能源 和 岩石 工程 .水利 岩石 工程 等 
宕 石 工程 重大 发 展 方向 的 需要 ,对 于 推动 DDA 和 
NMM 在 实际 问题 中 的 更 加 广泛 应 用 至 关 重 要 ,而 目 
前 这 方面 开展 的 工作 还 非常 有 限 。 建 议 在 DDA 和 
NMM 中 开发 完善 的 力 -水 - 热 多 场 耦合 功能 ,并 辅 以 
更 加 丰富 的 岩石 材料 及 结构 面 模型 以 及 岩石 破裂 模 
型 等 ,以 满足 岩石 工程 实际 问题 模拟 的 需要 。 

4) 三 维 方法 问题 

相 较 于 二 维 DDA 与 NMM 的 大 量 研究 与 成 功 应 
用 ,受到 接触 检测 与 处 理 、 计 算 效率 以 及 三 维 建 模 手 
段 等 因素 的 制约 ,三 维 DDA 和 NMM 的 发 展 与 应 用 


岩石 力学 与 工程 领域 广 受 认可 的 主流 数值 模拟 方 
法 。 本 研究 对 DDA 和 NMM 自 提 出 以 来 在 理论 与 算 
法 以 及 在 岩石 力学 与 工程 领域 的 应 用 研究 方面 的 工 
作 进行 了 较为 全 面 的 梳理 和 总 结 ,在 此 基础 上 总 结 
出 了 DDA 和 NMM 发 展 所 面临 的 主要 挑战 ,包括 计 
算 控 制 参数 .计算 效率 、 多 物理 场 、 三 维 方法 程序 的 
商业 化 开发 等 5 个 方面 ,并 给 出 了 应 对 建议 供 同行 
参考 ,以 期 DDA 和 NMM 的 更 好 发 展 与 应 用 ,更 好 服 
务 于 国家 岩石 力学 与 工程 模拟 重大 需求 。 
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